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Abstract 

We address the problem of evaluating an //-function when only a small 
number of its Dirichlet coefficients are known. We use the approximate 
functional equation in a new way and find that is possible to evaluate the 
i-function more precisely than one would expect from the standard ap- 
proach. The method, however, requires considerably more computational 
effort to achieve a given accuracy than would be needed if more Dirichlet 
coefficients were available. 

1 Introduction 

L-functions are central to much of contemporary number theory. Two celebrated 
conjectures, the Riemann Hypothesis and the Birch and Swinnerton-Dyer Con- 
jecture are about values of L-functions and were discovered as a result of the ex- 
plicit computation of the Riemann zeta function and the Hasse-Weil L-function 
associated to an elliptic curve, respectively. These L-functions are, respectively, 
of degree one and degree two and it is interesting to verify analogous conjectures 
about special values and zeros for higher degree L-functions. Conjectures such 
as Bocherer's Conjecture P3 [TS] and the Bloch-Kato Conjecture [3] are about 
the central values of degree four and degree three L-functions, and the Grand 
Riemann Hypothesis asserts that all nontrivial zeros of a "nice" L-function, of 
any degree, lie along the critical line. 

In addition to these conjectures, there are many other conjectures for the 
statistical behavior of L-functions, arising from random matrix theory and num- 
ber theoretic heuristics [101 1?1 [5] . One of the main reasons those conjectures are 
believable is that large-scale calculations of the value distribution and the zeros 
of L-functions yield data that support those conjectures. 
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The L-functions we consider here are associated to Siegel modular forms. 
Our examples will use the first non-lift Siegel modular form on Sp(4,Z). The 
form has weight 20 and is usually denoted Y20 Background information about 
Siegel modular forms can be found in (5J [TTJ |T5]. For this paper the relevant 
information is that a Siegel modular form is acted on by Hecke operators T(n), 
which have eigenvalues X(n). It is the eigenvalues A(p) and X{p 2 ) for p prime, 
which are used to define the L-functions associated to the modular form. For 
T20 the eigenvalues X(p) have been computed for p < 997, and the eigenvalues 
X(p 2 ) for p < 79 This data is available at [IT] . 

Associated to a Siegel modular form F is a sequence of L- functions, of de- 
grees 4, 5, 10, 14, 16, etc. There is an L-function L(s,F,p) of degree n for each 
n dimensional representation p of the dual group of PGSp(4), namely Sp(4,C). 
The degree 4, 5, and 10 L-functions are called, respectively, the spinor, stan- 
dard, and adjoint, and are denoted L(s, F, spin), L(s, F, stan), and L(s, F, adj). 



Proposition 2.1 taken from [5], summarizes the properties of those L- functions 
for a weight k Siegel modular form on Sp(4, Z). 

The degree 4 and 5 L- functions were shown by Andrianov pQ and Bocherer [3] 
to have an analytic continuation and satisfy a functional equation. The degree 
10 L-function was only recently shown by Pitale, Saha, and Schmidt [13] to have 
an analytic continuation and satisfy a functional equation. Those properties for 
the L-functions of degree 14 and above are still conjectural. 

1.1 Evaluating L-functions 

We are concerned with numerically evaluating L-functions. The standard tool, 
which is used in popular computational packages |141 [5] is the approximate 



functional equation. See Proposition 3.1 for the precise formulation. 



There are two main difficulties in evaluating high degree L-functions. The 
first is that if the L-function L(s) has degree d, evaluating L(~ + it) using the 
approximate functional equation requires ^s> (1 + |<|) d ^ 4 Dirichlet series coeffi- 
cients. Here the implied constant depends on the L-function and the desired 
precision in the answer. For example, estimating the implied constant from the 
calculations in Section |4j to find the first 1 million zeros of L(s, Y20, adj), the 
degree 10 L-function associated to T20, would require using the approximate 
functional equation with around 1 trillion Dirichlet series coefficients. Note that 
1 million zeros is not even a large sample; for example it is probably not sufficient 
for testing various conjectures about the lower-order terms in the distribution 
of zeros. 

The second difficulty is that current methods are incapable of producing a 
large number of Dirichlet coefficients of the standard and adjoint L-functions 
of a Siegel modular form. The Fourier coefficients of index up to 3000000 have 
been computed for T 2 o H2 ■ These Fourier coefficients are used to compute the 
Hecke eigenvalues. Examination of formulas on pages 387 of [T^ shows that to 
find the eigenvalue A(n) of T(n), for n — p 2 , requires the Fourier coefficients of 
index up to n 2 — p 4 . 
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It gets worse. By (2.2 1 and ( 2.4 1 , the pth Dirichlet coefficient of the standard 
or adjoint L-function requires both X(p) and A(p 2 ). Thus, to determine the first 
n Dirichlet coefficients of those L-functions requires Fourier coefficients of the 
Siegel modular form of index up to approximately n 4 . The extensive calculations 
in are not even sufficient to determine the 83rd Dirichlet coefficient of the 
standard or adjoint L-functions of Y20 

Of course, one could determine more Dirichlet coefficients by first finding 
more Fourier coefficients of the cusp form. But the n A relationship makes this 
quite expensive, so with current methods it is not feasible to determine many 
more Dirichlet coefficients than currently known. It is possible that new meth- 
ods will be developed to determine the Hecke eigenvalues without extensive 
computation. The real problem will still remain: how to compute high-degree 
L-functions without requiring an enormous number of Dirichlet coefficients. 

That brings us to the theme of this paper: how accurately can one compute 
an L-function given a limited number of Dirichlet coefficients. As the above 
discussion indicates, this is a practical problem and there are many L-functions 
for which it is not currently possible to determine a reasonably large number of 
Dirichlet coefficients. 

As we describe in this paper, even without knowing many Dirichlet coeffi- 
cients, we were able to evaluate the L-functions to surprisingly high accuracy 
(surprising to us, anyway). In fact, we were not able to establish that there 
is an absolute limit to the accuracy one can obtain from only a limited num- 
ber of coefficients. However, our method is computationally expensive - much 
more expensive than evaluating the L-function in a straightforward way if more 
coefficients were available. If one could find an efficient way to determine the 
unknown parameters in our method, that could make it possible to quickly 
evaluate high-degree L-functions. See Section [6] for a discussion. 

In the next section we describe the L-functions we consider here, and in 
Section [3] we recall the approximate functional equation and how it is used 
to compute an L-function. In Section [3] we state the underlying problem and 
then describe our experiments evaluating L(s, T 2 o, stan) and L(s, T 2 o, adj). In 
Section [5] we describe a second method that performs a little better than the 
method in Section [4] and show how this second method can be used to approxi- 
mate unknown Dirichlet series coefficients. 



2 The L-functions 

The L-functions associated to a Siegel modular form are most conveniently 
described as Euler products. The local factors of the Euler product can be 
expressed in terms of the Hecke eigenvalues \(p) and A(p 2 ), but it is more 
convenient to express them in terms of the Satake parameters, ao,pt a i,pi an< i 
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a 2 . p , given by 



P 



,2fc-3 



a 2 aia 2 

a^a\a 2 + a^aia 2 + a^a\ + a^a 2 

2 2 2 22 22 2 

a a 1 Q! 2 + aoC^ + a a 2 + a , 



(2.1) 



A 



B 



where 



\(p) 2 = 4p 2k - a + 2A + B 

\(p 2 ) = (2-l/p)p 2k - 3 + A + B. 



(2.2) 



We suppress the p on the Satake parameters when clear from context. 

We rescale the Satake parameters to have the so-called "analytic" normaliza- 
tion = 1, a§aia2 = 1, which is possible if we assume the Ramanujan bound 
on the Hecke eigenvalues. This corresponds to a simple change of variables in 
the -L-functions, so that all our L-functions satisfy a functional equation in the 
standard form sHl-s. 

As an error check for the reader who may wish to extend our calculations, for 
T 20 we have A(2) = -840960, A(4) = 248256200704, and the Satake parameters 
at 2 are approximately 

{a ,ai,a 2 } = 



{-0.901413 + 0.43296i, 0.630904 - 0.775861i, -0.211226 + 0.977437i}- 



Here and throughout this paper, decimal values are truncations of the true 



Proposition 2.1 Suppose F e M fe (Sp(4,Z)) is a Hecke eigenform. Let qzq p, 
a i,p, &2, P be the Satake parameters of F for the prime p, where we suppress the 
dependence on p in the formulas below. For p € {spin, stan, adj} we have the 
L-functions L(s, F, p) := n lme Q p (p' s ,F, p)" 1 where 



Qp{X, F, spin) := (1 — otQX)(l — aootiX)(l — ctQa 2 X)(l — a^aia 2 X), 
Q P {X, F, stan) := (1 - X){1 - ai X){\ - a^X)(l - a 2 X){l - a^X), 
Q p (X,F,adj) := (1 - X) 2 (l - ai X){l - a^X)(l - a 2 X)(l - a^X) 



give the L-series of respectively, the spinor, standard, and adjoint L-functions. 



(2.3) 



values. 



(1 - ai a 2 X)(l - a^ 1 a 2 X)(l - a^X)^ - a^a^X), 

(2.4) 
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These L-functions satisfy the functional equations: 

A(s, F, spin) := T c (s + |)r c (s + k- §)L(s, F, spin) 

= (-l) fc A(l- S ,F,spin), 
A(s, L, stan) := r R (s)L c (s + fc - 2)r c (s + k- l)L(s, F, stan) 
= A(l - s,F, stan), 
A(s, F, adj) := r R (s + l) 2 r c (s + 1) 

x r c (s + fc - 2)r c (s + fc - l)r c (a + 2fc - 3)L(s, F, adj) 
= A(l-s,F,adj). (2.5) 



In (2.51, we use the normalized T-functions 

r R («) := ^- s/2 r(s) and T c (s) := 2(27r)- s r(s). 

The degree of an L-function is r + 2c where r and c are the number of Tr and Tc 
factors in the functional equation, respectively. The spin, standard, and adjoint 
L-functions described above are of degree 4, 5, and 10. The Ramanujan bound 
for a degree d L-function with Dirichlet series X) n >i ^« n_s i s given by: 

M<( d+ j _1 )- (2-6) 

Note that this is equivalent to the assertion that the Satake parameters sat- 
isfy |0!j,p| < 1. 

3 The approximate functional equation 

In this section we describe the approximate functional equation, which is the 
primary tool used to evaluate L-functions. The approximate functional equation 
involves a test function which can be chosen with some freedom. This will play 
a key role in our calculations. 

3.1 Smoothed approximate functional equations 

The material in this section is taken directly from Section 3.2 of 1 14] . 
Let 

Hs) = E £ (3-D 

n=l 

be a Dirichlet series that converges absolutely in a half plane, 5R(s) > <j\. 
Let 

A(a) = Q s (fln^s + Xj) j L(s), (3.2) 
with Q, Kj £ R + , K(Aj) > 0, and assume that: 
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1. A(s) has a meromorphic continuation to all of C with simple poles at 
8i,...,sg and corresponding residues r\, . . . , re- 



2. A(s) = eA(1 - s) for some e G C, |e| = 1. 

3. For any <t 2 < 03, L(cr + it) = (^(expt" 4 ) for some A > 0, as |t| — > 00, 
C2 < cr < 03, with ^4 and the constant in the 'Oh' notation depending on 
02 and 03. 



Note that (3.2 1 expresses the functional equation in more general terms 
than (2.5 1, but it is a simple matter to unfold the definition of Tr and Y^. 

To obtain a smoothed approximate functional equation with desirable prop- 
erties, Rubinstein [2] introduces an auxiliary function. Let g : C —> C be an 
entire function that, for fixed s, satisfies 

|A(> + s)g(z + s)z- 1 \ -> (3.3) 

as \^sz\ — » 00, in vertical strips, — xq < 3?z < xq. The smoothed approximate 
functional equation has the following form. 

Theorem 3.1 For s ^ {si, . . . , si}, and L(s), g(s) as above, 
A( s )g(s) = £ ^™ + s E ^/i ( s > ») + zQ 1 -* E^/2(l- S , ») (3.4) 



s — 5^ ^ — ' n* * — ' n x 

k—l n—1 n—1 



/!(«,„):=_ / J]r(K 3 (z + .s) + A 3 )z- 1 3 ( S + z)(Q/n) z dz 

-1 rv+ioo a 

/ 2 (l- s ,n):=— / J] r ( K ^ z + 1 - s )+^) ;z " 1 5( s - z )(Q/«)"^ 

(3.5) 

f > max{0, — 3?(Ai/ki + s), . . . , — 5ft(A a /K a + s)}. 



In our examples L(s) continues to an entire function, so the first sum in (3.4 1 
does not appear. For fixed Q,n, X,e, and sequence b n , and g(s) as described 
below, the right side of (3.4) can be evaluated to high precision. 

A reasonable choice for the weight function is g(s) = e tbs+cs , which by 
Stirling's formula satisfies (3.3 1 if c > 0, or if c = and \b\ < ird/4, where d is 
the degree of the L- function. Rubinstein |14| uses such a weight function with b 
chosen to balance the size of the terms in the approximate functional equation, 
minimizing the loss in precision in the calculation. In this paper we exploit 
the fact that there are many choices of weight function, and so there are many 
ways to evaluate the L-function. We combine those calculations to extract as 
much information as possible from the known Dirichlet coefficients. This idea 
is described in the next section. 
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In the computations we carry out below, we find it more convenient to use 
the Hardy Z-function in our computations instead of the L-function itself. The 
function Z associated to an L- function L is defined by the properties: Z(^+it) 
is a smooth function which is real if t is real, and \Z(\ + it)\ = + it)\. 

4 Exploiting the test function in the approximate 
functional equation 



If we let g(s) = 1 and s = \ + lOi in the approximate functional equation (3.4 1 
for the standard (degree 5) L- function of T20, we get 

Z{\ + 10i, T 20 , stan) = - 1835.424 - 395.011 b 2 + 1012.179 b 3 + 1906.603 64+ 
+ 2226.503 65 + ■ • • + 6.840 x 10~ 9 b 82 + 
+ 5.132 x 10~ 9 b S3 + ■ ■ ■ + 3.205 x 10~ 16 6 149 
+ 2.564 x 10" 16 6i 50 + --- . (4.1) 

If instead we let g(s) — e~ 3ls / 2 and keep s = \ + lOi then we have 

Z(\ + m, T 20 , stan) = 1.66549 + 1.39643 b 2 - 0.658439 63 + 0.726149 64+ 
- 0.88227 65 + • • ■ + 1.532 x 10~ 8 6 82 
+ 1.271 x 10~ 8 b S3 + ■ ■ ■ + 3.514 x 10~ 14 6149 
+ 2.309 x 10" 14 6150 + ••• ■ (4.2) 



Note that neither of the above expressions appears optimal: the first involves 
large coefficients, which will lead to a loss of precision. In the second, the terms 
do not decrease as quickly, so one must use more coefficients to achieve a given 
accuracy. An observation that we exploit is the fact that the above are just two 
among a large number of expressions for the value of the L- function at | + lOi. 

Recall that by |12l ITT] we know the Satake parameters of T 2 o for all p < 
79. Therefore we recognize two types of terms in the approximate functional 
equation, as illustrated in the above examples. There are terms for which we 
know the Dirichlet coefficients, such as 65, b$ 2 , and 6150- And there are terms 
with an unknown Dirichlet coefficient, such as 6s3 or &149- Actually, there is a 
third type of term, such as 6 166 — b 2 b S3 which is "partially unknown". We can 
estimate the unknown terms by applying the Ramanujan bound to the Dirichlet 
coefficient, and evaluate everything else precisely. Thus, once we choose a test 
function, we can evaluate an i-function at a given point as 

Z(s) = calculated_ value (s) ± error_estimate(s), (4-3) 

where both the calculated value and the error estimate are functions of the test 
function and the set of known Dirichlet coefficients. For later use, we write 

error_estimate(s) = S n (g,s)b n (4.4) 

n: b n unknown 
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where the S n (g, s) is the coefficient of b n in (3.4|. The product of S n (g, s) and 
the Ramanujan bound of the L- function for the coefficient is an upper bound 
for the error contributed to computation by the unknown coefficient b n . 



get 



Using the known b n and applying the Ramanujan bound (2.6 1 to (4.1 1 we 
Z{\ + KK, T 20, stan) = 3.0393070838088 ± 3.12 x 10~ 8 . (4.5) 



And for (4.2) we get 

Z{\ + lOi, T 20 , stan) = 3.0393070808076 ± 7.10 x 10" 



(4.- 



In Figure [l] we show the calculated value and error estimate for 



lOi, T20, stan) when evaluated with test functions of the form g(s) 
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Figure 1: The solid line is the calculated value and the dashed line is the error 
estimate in computing Z(^+10i, T 2 o, stan) using the available Dirichlet coefficients 
with the weight function g(s) = e~ l/3s where /3 is given along the horizontal axis. 
The vertical axis is log 10 of (the absolute value of) the actual value. For quite a 
wide range of test functions, the value of the Z-function at h + lOi is determined 
with some accuracy, achieving around 10 decimal digits of accuracy with the optimal 
choice of /3. 



As Figure [T] shows, there is a wide range of j3 for which it is possible to 
determine A(| + lOi, T 2 o, stan) with some accuracy. With the optimal choice 
of /? the error estimate is approximately 4 x 10~ 10 . 

Figure |2]shows the calculated value and error estimate for the adjoint (degree 
10) i-function Z(^ + 5i,T2o,adj) when evaluated with test functions of the 
form g(s) = e -^+mo(^)\ 

As Figure [2] shows, every test function of the given form leads to an error 
which is larger than the calculated value. Thus, we can determine that \Z(h + 
5?,T2o,adj)| < 0.25, but with individual test functions of the given form we 
cannot even determine if the actual value is positive or negative. 
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Figure 2: The analogue of Figure[l]for Z(| +5i, T 20 , adj) with the weight function 
g(s) — e~ %l3s+ ~^° ( s_51 ) where f3 is given along the horizontal axis. 



We now introduce a new idea for increasing the accuracy of these calcula- 
tions. 

4.1 Optimizing the test function 

In Figure [T] we see that there are many values of the parameters in the test 
function which give reasonable results. If there is some degree of independence 
in the errors, then there is hope for obtaining a smaller error by combining 
the results of those separate calculations. Write Z(s, T 20 , stan, /?) for the out- 
put of the approximate functional equation with weight function g(s) = e~ l @ s . 
Consider 

j 

Z(s,T 20 ,stan) = J] c fe Z(s, Y 2 p, stan, /3j) (4.7) 

where ^ cp j = 1. We make the specific choices 
s = | + lOi 

(Pi, Pi, Pa, Pi, Pb) = (lOJ 5' 10' 5' 2) 

(c 01 ,cp 2 ,c 03 ,c^,c 05 ) = (0.03150,0.18061,0.36563,0.31421,0.10801). (4.8) 

Recall that all decimal numbers are truncations of the actual values; one re- 
quires much higher precision than the displayed numbers in order to obtain the 



answers below. With the choice in (4.8 1, after substituting the known Dirichlet 
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coefficients and then using the Ramanujan bound, we find 



Z(s, T 2 o,stan) 



3.0393070864895278277801 + 2.688 • 10~ 19 b 83 + 1.147 • 10~ 16 6i 07 - 



+ 5.291 ■ 10" 18 6i37 + ■ • • + 1.216 • 10" 23 &i99 

-1-15 



=3.039307086489527827 ± 4.737 • lO" 15 . (4.9) 

Thus, by averaging only 5 evaluations of the L-function, the error decreased by 
a factor of almost 10~ 5 . 

The weights cp. in (4.8 1 were determined by finding the least-squares fit to 

Ram(b n ) 2 (y2c0A( 3 3A + m )) =0. ( 4 - 10 ) 

n: b„ unknown 7 



where Ram(b n ) is the Ramanujan bound (2.6 1 for b n , subject to X) c ft = 
Note that in our actual examples the vast majority of unknown coefficients 
have prime index, so the Ram(b n ) weighting is not important, but we include 
it for completeness. 

The error estimate in (|4.9|) is an L 1 estimate, not the L 2 estimate as shown 



in (4.10), so actually it is possible to choose slightly better weights that those 
used in our example. In Section[5]we show how to obtain the optimal result that 
can arise from combining different evaluations of the L-function, but for now 
we merely wish to illustrate the seemingly surprising fact that appropriately 
combining several evaluations can vastly decrease the error. 



4.2 Results 

In Figure [3] we plot the error obtained by combining varying numbers of weight 
functions, where we evaluate Z(^ + lOi, T20, stan) with a weight function g(s) = 
e l/3s with j3 = j/10 with —10 < j < 25. The horizontal axis is the number of 
terms averaged, where we start with /3 = | and first use those /3 which are 
closest to i. The vertical axis is the error estimate on a log 10 scale. The lowest 
point on the graph, when we average all 36 evaluations, corresponds to 

Z(i+10i,T 20 ,stan) = 

3.0393070864895284810824603284422250910 ± 2.79 x KT 35 . 

It is not clear from Figure [3] whether one would expect to obtain an arbitrar- 
ily small error by combining sufficiently many test functions in the approximate 
functional equations. This is discussed in Section [6] 

We briefly describe the corresponding calculations for Z(h + 5z,T2o,adj). 
Recall that, as shown in Figure [2] with a single test function of the standard 
form we are not able to determine whether that value is positive or negative. 
Now we combine five evaluations in the analogous way: 

.7 

Z( s ,T 20 ,adj) =^ % Z(s,T 20 ,adj,/3 J -) (4.12) 
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Figure 3: The error obtained from a least-squares minimization of the error for 
combining n evaluations of Z(~ + lOi, T 2 o, stan) using the weight functions g(s) = 
e~ ll3s with /? = j/10. The horizontal axis is n and the vertical axis is log 10 of the 
resulting error estimate. 



where ^ cp. = 1. Here the weight function is g(s) = e % P s +i s 5l ) / 50 °. 
We make the specific choices 

s = \ + 5i 

(ft, ft, (f.f^f.f' 3 ) 
{c p ^c h ,cp 3 ,c^,cp B ) = (0.035863,0.33504,0.47934,0.13827,0.01146). (4.13) 

The result is 

Z{\ + hi, T 20 , adj) = 0.01556 ± 0.0049. (4.14) 

So we have determined that Z(h + 5i,T2o,adj) is positive, but we can only be 
certain of one significant figure in its decimal expansion. Using this method, 
the best result we were able to obtain, by averaging 11 evaluations, is 

Z{\ + hi, T 20 , adj) = 0.01558768 ± 0.00016. (4.15) 

Averaging more weight function actually makes the result worse. The expla- 
nation is simple: since the adjoint L-function has high degree, the error terms 
5 n decrease very slowly. The least-squares fit does not properly take into ac- 
count the contribution of a large number of small terms, so the least-squares fit 
actually has a large L 1 norm. 

In Section [5] we introduce a different method which avoids these shortcom- 
ings. 
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5 Linear programming 



In this section we view the evaluation of the L-function as an optimization prob- 



lem. For example, we can view the equality of the expressions in (4.1 1 and (4.2 1 
as a constraint on the value of the L-function. Thus, the same calculations 
which were used as input for the least-squares method described in the previous 
section can also be used as input to a linear programming problem. 

We set up the linear programming problem as follows. Let Z(^+10i, T 2 o, stan, gj) 
denote the evaluation of Z(i + lOi, T 2 o, stan) using the weight function gj in 
the approximate functional equation. One evaluation, Z(| + lOi, T 2 o, stan, gi), 
is taken as the objective. The other evaluations are taken in pairs and 

Z{\ + 10i,T 20 , standi) - Z{\ + lOi, T 20 , stan, gj) = (5.1) 

is interpreted as a constraint. The other constraints come from the Ramanujan 
bound on the unknown coefficients. 

In the above description there are infinitely many unknowns and constraints. 
We eliminate the unknown coefficients b n with n > 1000 by using the Ramanujan 



bound, replacing the equality (5.1 1 by a pair of inequalities. Thus, we have a 
straightforward linear programming problem, which we use to determine the 
minimum and maximum possible values of Z(~ + lOi, T 2 o, stan). 

We implemented this idea using the same set of test functions we used for 
the least-squares method. The calculations were done in Mathematica, using 
the built-in LinearProgramming function with the Method -> Simplex option. 
In all cases the linear programming method gave better results, but not spec- 
tacularly better. Figure [4] shows the ratio of the errors in the results of the two 
methods, on a log 10 scale. For example, when using 30 equations the error from 
the linear programming approach was approximately 1/10 the error from the 
least-squares method. 



5.1 Computing unknown coefficients 

An interesting side-effect of the linear programming approach is that it allows 
us to obtain information about the unknown coefficients. Instead of treating the 
value of the L-function as the objective, we can use an unknown coefficient as 
the objective. Note that all the other constraints in the problem are unchanged. 
Using the same method as described above we find the following coefficients of 
the standard L-function of T 2 o: 

b S3 = 0.4884558312724 ± 2.4 x 10~ 12 
6 89 = 0.10561760640 ± 2.7 x 10~ 10 

697 = 0.468135808 ± 1.5 x 10~ 7 (5.2) 

These results can be checked once more Fourier coefficients of T 2 o are com- 
puted. By the n 4 argument given in in the Introduction, computing exact values 
of bp for p < 97 will take about twice as much work as it took to compute those 
for p < 79. 
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Figure 4: The ratio of the errors in the linear programming and least-squares meth- 
ods for evaluating Z{\ + lOi, T20, stan) using n equations. The horizontal axis is 
n and the vertical axis is log 10 of the ratio of errors. 



6 Conclusions and further questions 

We have shown that, at the cost of a lot of computation, one can evaluate an L- 
function to high precision using only a small number of coefficients. That this is 
theoretically possible is not surprising: L- functions are very special objects, and 
the data we have for the L-function considered here (the functional equation and 
the first several coefficients) presumably specify the L-function uniquely. Thus, 
in an abstract sense there is no new information in the missing coefficients. But 
the question remains as to whether our methods accomplish this in practice. 

Question 6.1 Can the method of calculating an L-function by evaluating the 
approximate functional equation and then averaging to minimize the contribu- 
tions of the unknown coefficients, determine numerical values of the L-function 
to arbitrary accuracy? 

Because the approximate functional equation requires a huge number of 
terms to evaluate a high-degree L-function, it would be significant if the weights 
we obtained by our methods could be determined without actually calculating 
all the terms with unknown coefficients. 

Problem 6.2 Devise a method of determining an optimal weight function in 
the approximate functional equation without first calculating a large number of 
terms which do not actually contribute to the final answer. 

Implicit in the above problem is the requirement that one knows the calcu- 
lated value as well as an estimate of the error. 

Problem 6.3 Is there any meaning to the weights determined by the least- 
squares method? 
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The weights which appear in our least-squares method depend on the point 
at which the L-function is evaluated. It might be helpful to consider a case 
where there are two L-functions with the same functional equation, such as the 
spin L-functions of Sp(4, Z) Siegel modular forms of weight k > 22. 

Without progress on these problems, or a completely new method, there is 
little hope of making extensive computations of high-degree L-functions. 
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